Adventures in Classical PSF 
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Need for Coronagraphy 


• Circumstellar Disks, exoplanets, stellar 
companions are often inconveniently dose to a 
bright object (host star) 

• Exposing sufficiently deeply to detect the object 
of interest can mean that you overexpose the 
instrument you are using and swamp the signal of 
interest 

• his talk will focus on using a simple coronagraph, 
and how best to separate the signal of interest 
from light from the star 
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The PSF 


Optical systems, and HST is no exception, 
typically spread the light from an unresolved 
source due to diffraction, scattering in the 
telescope, and in the science instrument, and 
in some cases within the detector system. 

"his is termed the point spread function PSF). 

For the majority of circumstellar disks and 
exoplanets signal in the wings of the 
PSF»signal of interest. 


Classical PSF Subtraction 


• Simplest of the techniques that will be 
covered in this hands-on demonstration 

• Use a suitably chosen other observation as an 
estimate of the light from the star that you 
want to get rid of. 

• Need to match the science observation in 
terms of factors affecting the shape of the PSF, 
and those affecting temporal variation in the 
measured PSF. 


Shape of HS" PSF depends on Source 
Spectral Energy Distribution 


Targets for coronagraphic observation are 
typically on the bright side for direct 
Imaging. 

STIS CCD used for coronagraphic imaging 
Has throughput from 0.2-1. 0 microns 

Effective wavelength of the image is 
a strong function of Teff. 

Here we see the observed response for 
a white dwarf (Feige 110) - bulk of 
signal is at wavelengths <5000 A 
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A1V simulation 


By A1V, the bulk of the signal is in the 
Conventional optical, but with broad 
Wings. 
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M2 Simulation 


By early M, the effective wavelength 
of the STIS image is 7700 A. 

Sensitivity to color of source decreases 
as bandwidth decreases, but is typical 
of all 3 HST coronagraphs. 

JWST will have other issues, namely 
a sensitivity to thermal emission from 
The inner disk (IR excess), and the fact 
that diskless stars will be blue 
compared to any system with an IRE. 



Sagan Summer Wo rks h o o 2 0 14 


And also on "actors affecting focus 


• HS is in a low-Earth orbit and experiences 
changing thermal conditions. 

•Scheduling for HST does not include thermal 
effects 

•differences between focus conditions for the 
target and what you are using as an estimate of 
the PSF result in radial streamers - differences in 
the dispersed speckles 
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:ase 1: Use the same star, different roll 

of spacecraft on sky 


AU Mic - obiw36030jflt.fits Subtracting obiw35030_flt.fits 
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Introduction to IDP3 


• For exploring PSF subtraction you need a display-oriented 
tool that will let you play with registering, scaling, and 
subtracting one image from another, while allowing you to 
tweak the display scaling, color, etc. to your preferences. 

• Various tools used by the HST IDTs, and other teams, but 
we have chosen IDP3 developed by the NICMOS team 

• Starting point is file in the detector frame with NO 
geometric distortion correction or mapping to the sky. 

• Get into IDL, type IDP3 

• Large window should appear with File, Images, Adjust, Edit, 
Spitzer, and Help on menu bar 
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Loading Data 


Choices 

• STIS data are typically multiple 
sub-exposures to facilitate 
cosmic ray removal. 

• Can load the first image of a 
set, or the full set. 

• Under file, go to the load 
images item 

• To read the first readout, go to 
load image on the 2 nd menu, 
to read in a set, go to Load 
Multiaccum (NICMOS 
terminology). 
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Selecting the file to be loaded 


hi Iter 



Directories 



+ * 




hiles 

*alignwcs + txt 

| aum i c_36m35 * f i ts 
hdl81327_co 1 ormatchN I C * png 
hdl81327_mosa i c_6 * png 
hdl81327_r i ng_st is*png 
\ hdl01327ul2040mhdl34970 + png 
\ hdl81327yl2040mhr4413 * png 
jj hdi81327vl2040mvarpsf s + ai 
l idp3_tmp + txt 


ob i hit 1040 _f 1 t 4 . f i t-5 


| ob i wl3040_f 1 1 + f i ts 
I obiwl7040_f lt*f its 
I ob i wl8040_f 1 1 + f i ts 
ob i w35030_f 1 1 * f i ts 
3 ob i w38030_f 1 1 + f I ts 
I ob i w93040_f 1 1 + f i ts 
! ob i w97040_f 1 1 + f i ts 


Select on 
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Loading multiple reads 


File>load mage>load 
multiaccum 

v Please SeJeck Image “ £ 

sctory 

Dlunes/datal/Users/cgrady/Sagan.data.Grady/f 


To load all reads from an image, type *; or you 
can specify individual reads by number 
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Files 

alignucs.t^t 
auniic_36ra35,f its 
hdl81327_colooiatchNIC, png 
hdl91327_mosaic_G .png 
hd!81327_r ing_stis *png 


hdl81327wl204U»hdl3497A +P ng 
hd!81327wl2w4nmhr4413*png 
hdl81327wl2u4i iwvar-psf s,a i 
idp3_tnp.txt 


ob iwiiW.rsedl an+f its 
obiwl3n4"_f lt*f its 
obiwl704^_flt,fits 

obiw3S03u_f It .fits 
obiwSGOSAjflt.fits 

obiu93Q40_flt„fits 

t 


Select reads to load,.... 

y 

jj Syntax: * = all # = single 


action 

iulltJ40_flt.fit& 
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Showing the images 
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Adjusting the Display 


• You can adjust the image dynamic range, color 
table, and scaling (linear, log, square root) 
using adjust> display menu item. 

• HST coronagraphic images typically have a 
large dynamic range, so I find log scaling 
works well. 

• You can also shrink the window to the actual 
data size - IDL increments from the lower left 
image corner- using the resize display option. 
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Image registration 
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image Registration 


• In show images window select the image to be 
subtracted as the one with * 

• Now go to the Adjust>adjust position menu 

• Default increment is 1.0 pixels 0.05 ’ for STIS), 

• I find 0.05 pixel works better 

• Use the multiple arrow panel to "drive" the image to 
be subtracted so that there is not a significant 
dark/light asymmetry in the image and the diffraction 
spikes are largely nulled 

• In this case since we were using data for the same star 
taken in consecutive orbits, we did not need to adjust 
the flux scaling for the net image. 
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Saving Your Image 
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Saving Your image 


• You can save for use within IDP3 

• Or write a file for external use or both 

• Go to File>save display menu, give a file name 
and path. 

• Congratulations you now know the basics of 
classical PSF subtraction. 
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Now for a face-on disk... 



3 aran Smi roer Vtorksh 



Strengths and Weaknesses of ADI-like 

strategies 

• Best for point source detection and edge-on 
disks such as AU Mic 

• Can eliminate any signal which is azimuthally- 
symmetric over the roll angle - face on disks 
can be eliminated. 

• Can end up with a mess if the nebulosity is 
very structured 

• Need alternate robust technique for removing 
the PSF, which conserves flux 
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Case 2: Non-Contemporaneous PSF 
template observations 

• Science target HD 18132 , mid-F member of (3 
Pic moving group 

• Disk marginally resolved by Herschel 
(Lebreton et al. 2012) 

• Disk resolved in scattered light by NICMOS 
and ACS (Schneider et al. 2006) 

• Location of bulk of debris constrains SED 
modeling; asymmetries can constrain planets. 
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Using archival data for HR 4413 as 

template 



HD 181327vl2-hr4413 
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sing HD 1349 0 as template 





Sommer Workshop 2014 



What you get rom such processing 


• Detection 

• Bulk of signal in ring - inner edge, outer edge, 
inclination, can compare with predictions for location 
of debris belt from FIR data assuming that the grains 
are large, compact grains (gray or blackbody). 

• Exterior to ring, additional, azimuthally asymmetric 
nebulosity 

• There are residuals which depend on choice of PSF 
template data, and this is a bright disk. Residuals 
become more of a nuisance as the surface brightness 
of the disk decreases. 
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ase 2 hands-on 


• Load files obiwll040_flt.fits 

• obiwl2040_flt.fits 

• obiwl3040_flt.fits - psf template 

• obiwl7040_flt.fits - psf 

• obiw93040_flt.fits - psf 

• obiw97040jflt.fits - psf 

• Now try subtracting 12040 from 11040 and using 
the PSF template data, some planned and some 
template used for another target in this program 
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Case 3: Color-matched, 
contemporaneous PSF data 
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Comments on color-matched, 
contemporary PSF subtraction 

• Suppress residuals, since have allowed HST to 
come to a quasi-equilibrium 

• Note large number of hot pixels, cosmic ray 
events, etc. - these can be removed by 
combining a suite of science target-PSF data 
where HST is rolled between observations - 
reduce S‘ r IS wedge to a quasi-circular occulted 
zone with r=0.35 , and can median filter to 
remove hot pixels. 
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HD 181327 - smaller obscuration 



Image Credit Schneider et al. 2014 
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HD 18132 after merging images 
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Deprojection and compensation for r 2 

illumination gradient 
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Further improvements & Summary 


• Can largely remove the remaining residua s 
using filtering techniques 

• Now are at the point that you can begin 
science 

• PSF subtraction with HST requires choice of 
suitable template targets, planning the 
observations so that the template is taken as 
close in time to the science data as feasible, 
and straightforward data reduction. 
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